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ABSTRACT 


A  method  is  described  for  determining  the  phase  constant  of  the  rock  me¬ 
dium  below  the  surface  of  the  earth.  The  constant  is  deduced  from  the  meas¬ 
ured  phase  angle  of  the  mutual  impedance  between  two  vertically  polarized 
antennas  placed  in  drill  holes  extending  into  the  rock. 

Experiments  were  conducted  on  Cape  Cod,  Massachusetts,  using  the  meth 
od  to  determine  the  phase  constant  of  the  rock  between  two  1000-foot  drill 
holes,  separated  by  a  distance  of  6000  feet,  at  frequencies  between  100  and 
4200  cps.  Electrically  short  monopoles  were  used  in  the  drill  holes,  with  the 
overburden  above  the  rock  serving  as  a  ground  plane.  Reference  phase  was 
derived  from  Loran  C  transmissions  at  both  the  transmitting  and  receiving 
ends  of  the  path. 

The  measured  phase  constants  at  these  frequencies  indicate  a  large  loss 

tangent  for  the  medium.  The  bulk  conductivity  was  thus  deduced  to  be  1.1  x 

-3 

10  mhos/meter,  and  the  measurements  did  not  yield  an  estimate  of  the  rela¬ 


tive  dielectric  constant. 


I.  INTRODUCTION 


The  rock  strata  below  the  surface  of  the  earth  is  a  dielectric  medium 
which,  although  lossy,  can  support  propagation  of  radio  waves.  *  The  medium 
is  generally  inaccessible  except  through  deep  drill  holes  in  which  antennas  for 
transmitting  and  receiving  can  be  inserted.  Since  the  electrical  characteristics 
of  the  rock  strata  vary  from  location  to  location,  the  complex  phase  constant 
k  =  p  -ja  of  any  transmission  path  between  two  antennas  must  be  estimated 
either  through  surface  measurement  or  by  measuring  the  transmission  loss 
of  the  path. 

The  transmission  loss  or  attenuation  method  has  been  successfully  used 
to  evaluate  the  attenuation  constant  a  from  which  the  conductivity  of  the  medi¬ 
um  can  be  deduced  if  the  loss  tangent  is  high.  *  An  equally  important  method 
of  evaluating  the  medium  characteristics  is  to  measure  the  phase  constant  p, 
i.  e.  ,  the  transmission  delay.  For  a  very  lossy  medium  the  transmission  de¬ 
lay  (or  phase  shift)  measurements  yield  only  the  conductivity.  When  the  loss 
tangent  is  small,  the  transmission  delay  is  a  measure  of  the  relative  dielec¬ 
tric  constant  of  the  medium;  whereas  the  transmission  loss  method  does  not 
result  in  an  independent  estimate  of  either  conductivity  or  the  relative  dielec¬ 
tric  constant. 

The  phase  method  described  in  this  report  is  based  on  the  measurement 
of  the  phase  angle  of  the  mutual  impedance  between  two  antennas  as  described 
in  Section  II.  A  practical,  simple,  and  convenient  way  of  providing  reference 
phase  through  the  use  of  Loran  C  transmission  is  shown  in  Section  III.  Sec¬ 
tion  IV  shows  the  experimental  results  of  a  test  conducted  on  Cape  Cod  in¬ 
volving  a  transmission  path  length  of  6000  feet  in  the  granitic  strata  below  a 
layer  of  overburden. 
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II.  THEORY 


In  conducting  deep  strata  propagation  investigations,  access  to  the  rock 
medium  is  usually  through  a  limited  number  of  deep  vertical  drill  holes.  Meas¬ 
urements  of  medium  characteristics  should  then  be  those  which  use  the  least 
number  of  access  holes.  For  in  situ  measurement  of  the  phase  constant  of  a 
propagation  medium,  two  drill  holes  are  required.  We  shall  show  that  the 
phase  constant  of  a  path  in  a  dissipative  medium  can  be  determined  from  the 
measurement  of  the  phase  angle  of  the  mutual  impedance  between  two  antennas 
placed  in  drill  holes  at  two  ends  of  a  path. 


If  a  linear  antenna  of  effective  length  h^  has  an  input  current  1^,  an  elec¬ 
tric  field  is  excited  in  its  surrounding  medium.  If  a  receiving  antenna  is  par¬ 
allel  to  the  transmitting  antenna  and  both  antennas  have  the  same  equatorial 
plane,  it  can  be  shown  that  for  an  effective  receiving  antenna  height  of  h^,  the 
mutual  impedance  between  the  antennas  is  given  by 


Z 

m 


hlh2W,\)  R  ,  ,  ,  1  ,  -jkR 

(jkR  +  1  +  j£p-)  e  J 


2  RZk 
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where 


co  =  radian  frequency 

|jlo  =  permeability  of  free  space  ^ 

k  =  iM  1  _j  jy) 

R  =  distance  between  antennas 

and  a  and  (5  are  respectively  the  attenuation  and  phase  constants  of  the  medium 
defined  as 


P  -ja  =  w 


-  J'P 


(3) 
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where 


£  =  relative  dielectric  constant 

r 


CO  £  £ 

o  r 


The  quantity  p  is  known  as  the  loss  tangent  and  c  is  the  conductivity.  MKS 
units  are  used  and  the  medium  is  assumed  nonmagnetic. 

If  we  write 


r-f  ry  ill 

Z  =  Z  e 
m  m 

the  phase  angle  0  of  the  mutual  impedance  may  be  obtained  by  substituting 
equation  (2)  in  equation  (1),  whence 


-PR  -h  2  tan"1  1  +  tan  _1  - - ^ 

I1  i 


-  p  -  (OR)  $2 


1  +  2  |5R 


0  is  a  function  of  u,  a  ,  and  e^.  We  observe  in  equations  (3)  and  (4)  that  p 
is  a  positive  real  quantity  so  that  the  ratio  -  is  limited  to  values  between  1  and 
0,  that  is 


The  ratio  will  assume  a  value  close  to  1  or  to  0  for  only  a  very  small  range 
of  frequency.  We  thus  may  simplify  equation  (5)  by  considering  pR  as  an  in¬ 
dependent  variable.  Figure  1  shows  the  relation  between  0^  and  pR  with  ^ 
as  a  parameter. 


It  is  seen  in  Figure  1  that  for  PR  greater  than  3  radians,  0^  is  essentially 
independent  of  Therefore,  if  0^  is  a  measured  quantity  less  than  -100 
the  phase  constant  of  a  path  may  be  determined  with  very  little  error  by  assum¬ 
ing  =  I,  F  or  lower  values  of  pR,  the  actual  ratio  p-  will  generally  be  a 
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function  of  frequency.  Furthermore,  since  the  slopes  of  the  0  vs  0R  curves 

m 

are  very  small  when  0^  is  small,  it  will  be  quite  difficult  to  determine  PR  with 
reasonable  accuracy.  Thus,  the  transmission  delay  method  is  primarily  use¬ 
ful  when  the  electrical  path  length  is  sufficiently  large  so  that  the  receiving 
antenna  is  in  the  far  field  zone  of  the  transmitting  antenna. 

It  should  be  noted  that  in  a  lossy  medium,  when  the  frequency  is  sufficient¬ 
ly  low,  p  will  be  nearly  unity,  so  that  0^  approaches  zero  as  PR  approaches 
zero.  This  is  an  important  feature  because  it  can  be  used  to  resolve  ambiguity 
of  an  integral  number  of  cycles,  which  may  not  be  otherwise  apparent  in  meas¬ 
uring  0  . 

6  m 

If  loss  tangent  p  is  small,  say  p  <  0.  6,  then 


p  *“Viv77 

(6) 

„  <rj  ^o 
°  “  2  V  EE 

(7) 

The  attenuation  constant  a  is  a  function  of  both  tr  and  e  and  is  independent  of 
frequency  while  P  is  only  a  function  of  e^  and  frequency.  Thus,  the  dielectric 
constant  of  the  medium  can  be  determined  only  if  both  attenuation  and  phase 
methods  are  used.  On  the  other  hand,  when  p  is  large,  say  p  >  10 


P  *  a  -  (8) 

The  phase  and  attenuation  constants  are  very  nearly  equal  ana  proportional  to 
the  square  root  of  tr  and  frequency,  but  independent  of  e  .  These  remarks  ap¬ 
ply  to  a  dissipative  medium  in  which  <r  and  e^  are  constant  with  frequency. 
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III.  INSTRUMENTATION 


Common  reference  phase  must  be  provided  at  receiving  and  transmitting 

ends  of  a  path  if  the  phase  constant  of  the  medium  is  to  be  determined.  In 

coastal  regions  of  the  United  States,  this  reference  is  readily  available  in  the 

form  of  Loran  C  transmission.  Loran  C  is  a  radio  system  primarily  intended 

2  3 

as  a  navigational  aid.  ’ 

Briefly,  Loran  C  signal  is  a  pulse-modulated  100-kc  carrier-frequency 
signal.  The  basic  repetition  rate  is  20  cps.  Modulation  and  carrier  frequen¬ 
cies  are  phase  coherent.  The  transmission  mode  of  Loran  C  for  short  ranges 
is  based  on  ground  wave  propagation.  The  rise  time  of  the  pulse  envelope  is 
30  microseconds,  so  that  by  observing  the  RF  cycles,  the  waveform  of  a  pulse 
may  be  used  as  an  accurate  time  marker.  Being  a  ground  wave  signal,  the 
received  waveform  of  Loran  C  is  delayed  from  the  transmitted  waveform  by 
an  amount  proportional  to  the  known  distance  between  the  receiving  site  and  the 
Loran  C  station. 

In  our  deep  strata  experiment,  two  receiving  sites  are  involved.  The  time 
base  at  each  site  is  derived  from  the  Loran  C  signal.  There  is  a  finite  dis¬ 
crepancy  between  the  time  bases  corresponding  to  the  difference  in  distances 
from  these  sites  to  the  Loran  C  station.  This  discrepancy  can  be  easily  ac¬ 
counted  for  if  one  so  wishes.  The  time  base  generated  by  the  Loran  C  signal 
is  relative,  in  the  form  of  time  markers  at  a  rate  of  20  pulses  per  second. 

The  block  diagram  in  Figure  2  shows  the  instrumentation  required.  A  100-kc 
receiver  is  used  to  amplify  the  Loran  C  signal.  The  output  of  the  receiver 
is  displayed  on  a  monitor  scope.  The  purpose  of  this  scope  is  to  generate  a 
more  convenient  waveform.  The  scope  is  internally  triggered,  and  the  trigger 
circuit  is  adjusted  so  that  the  sweep  of  the  scope  starts  at  the  desired  point  on 
the  waveform  of  the  Loran  C  pulses.  The  gate  output  of  the  scope  is  a  20-cps 
rectangular  wave  with  its  leading  edge  coincident  in  time  of  arrival  with  the 
phase  of  the  Loran  C  pulse.  If  this  same  instrumentation  is  used  at  both  test 
sites,  there  will  then  be  a  common  reference  phase  at  20  cps. 
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The  frequency  for  the  propagation  test  is  derived  from  the  20-cps  gate 
output  of  the  monitor  scope.  The  waveform  of  the  gate  output  is  rich  in  har¬ 
monics.  A  suitable  instrumentation  at  the  transmitter  end  of  the  deep  strata 
propagation  path  is  shown  in  Figure  3.  A  wave  analyzer  selects  the  desired 
harmonic  frequency  of  the  20-cps  gate  output  waveform.  The  selected  har¬ 
monic  is  available  at  the  output  of  the  wave  analyzer  (Hewlett-Packard  Type 
302A)  for  driving  the  transmitter.  A  current  transformer  in  series  with  the 
antenna  converts  the  current  waveform  into  a  voltage  waveform  which  is  ap¬ 
plied  to  the  vertical  input  of  a  measuring  scope.  The  measuring  scope  is  used 
to  determine  the  phase  of  the  antenna  current  relative  to  the  20-cps  reference 
phase.  For  this  measurement,  it  is  necessary  to  use  only  the  gate  output  from 
the  monitor  scope  to  trigger  a  delayed  sweep  on  the  measuring  scope.  The 
calibrated  delay  adjustment  on  the  scope  is  used  to  determine  the  phase  of  the 


At  the  receiving  end  of  the  path  a  measuring  scope  is  also  used,  as  shown 
in  Figure  4.  Here,  the  vertical  input  to  this  scope  corresponds  to  the  voltage 
induced  on  the  receiving  antenna.  In  this  case,  the  wave  analyzer  is  used  as 
a  narrowband  receiver.  The  phase  of  the  antenna  voltage  is  denoted  by  0^  and 
if  the  phase  shift  of  the  receiver  is  negligible,  the  measuring  scope  will  meas¬ 
ure  0^  directly.  Since  the  angles  0^  and  0^.  measured  by  the  scopes  are  in 
terms  of  phase  delay,  the  phase  angle  of  the  mutual  impedance  is  then 


0 

m 


(9) 


This  represents,  of  course,  an  idealized  situation.  In  practice,  corrections 
may  be  required  to  account  for  phase  shifts  in  transmission  lines  connecting 
the  antennas  and  for  phase  shift  in  the  receiver.  The  correction  can  be  readily 
measured  or  calculated  if  the  impedances  of  the  antennas  and  transmission 
lines  are  known. 
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IV.  CAPE  COD  EXPERIMENT 


The  instrumentation  described  above  for  phase  measurement  was  employed 
recently  in  tests  on  Cape  Cod,  Massachusetts.  The  drill  holes  are  located  on 
the  Town  Dump,  Brewster  and  Tubman  Road,  Brewster.  These  are  referred 
to  as  Brewster  and  Tubman  drill  holes,  respectively.  The  drill  holes  are 
6000  feet  apart  and  are  both  1000  feet  deep.  The  granitic  rock  in  that  area  is 
separated  from  the  surface  of  the  earth  by  a  layer  of  overburden  about  100  feet 
thick. 

Antennas  were  inserted  into  the  drill  hole  just  below  the  overburden,  which 
was  used  as  a  ground  plane  for  the  monopole  antennas.  The  antennas  used 
were  the  insulated  type  with  short-circuit  termination  and  low  impedances.  ^ 

The  test  frequencies  ranged  from  100  cps  to  4220  cps.  No  correction  for  phase 
shift  in  the  antenna  systems  was  found  necessary.  The  phase  shift  in  the  re¬ 
ceiver  amounted  to  90°  at  the  lowest  frequency  and  reduced  to  negligible 
amount  when  frequency  was  increased. 

The  measured  phase  angle  0^  (in  degrees),  determined  from  equation  (9), 
is  plotted  as  a  function  of  frequency  in  Figure  5.  It  is  observed  that  the  phase 
angle  0  can  he  extrapolated  to  the  origin  (0  =  0,  pR  s  0),  indicating  that 

there  is  no  ambiguity  of  multiples  of  360°  involved.  The  measured  phase  an¬ 
gles,  0^ ,  are  used  to  determine  (3R  from  the  curve  in  Figure  1  for  -j~  =  1.  The 
resultant  values  of  |3R,  in  radians,  are  plotted  in  Figure  6  as  a  function  of 
frequency.  A  straight  line  drawn  through  these  points  on  log-log  paper  has  a 
slope  of  0.  5,  indicating  that  (3  is  proportional  to  the  square  root  of  frequency 
as  in  equation  (6)  and  that  the  loss  tangent  is  large. 

One  may  then  deduce  from  these  measurements  that  the  average  conductiv- 

-  3 

ity  of  the  path  between  the  Tubman  and  Brewster  drill  holes  was  1.1x10 
mho/meter.  The  value  of  tr  deduced  from  the  attenuation  constant  was  very 
close  to  this  value.  Such  values  indicate  a  poor  medium  for  rock  propagation. 
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V.  DISCUSSION 


A  method  of  determining  the  real  phase  constant  (3  of  a  conducting  medium 
has  been  devised.  The  phase  method  is  useful  primarily  in  the  far  field  when 
the  phase  angle  of  the  mutual  impedance  between  two  antennas  is  essentially 
independent  of  the  loss  tangent  of  the  medium.  If  loss  tangent  is  large,  the 
phase  method  has  no  greater  utility  than  the  simpler  attenuation  constant  meth¬ 
od.  The  attenuation  method,  however,  must  be  implemented  by  the  phase 
method  to  yield  the  conductivity  and  relative  dielectric  constant  of  the  medium, 
if  the  loss  tangent  is  small. 

The  application  of  the  phase  method  is  particularly  simple  in  areas  where 
reception  of  Loran  C  is  good.  Otherwise,  more  cumbersome  means  must  be 
employed  to  provide  common  reference  phase  or  time  base. 

One  interesting  aspect  of  the  phase  method  is  that  the  transmission  is  phase 
coherent  with  the  reference,  so  that  arbitrarily  narrow  bandwidths  may  be  used 
in  detection  to  enhance  the  signal-to-noise  ratio.  This  not  only  increases  the 
usable  frequency  for  phase  measurement,  it  also  improves  the  reception  for 
attenuation  measurement. 
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Figure  1.  Phase  Angle  of  Mutual  Impedance 
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Figure  2.  Generation  of  20-CPS  Reference  Phase 
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Figure  3.  Block  Diagram  for  Phase  Measurement  at  Transmitting  Site 
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ure  4.  Block  Diagram  for  Phase  Measurement  at  Rec 
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